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ABSTRACT 

We report the discovery of a new hot-Jupiter, KOI-196b, transiting a solar-type star with an orbital period of 1.855558 days ± 0.6s 
thanks to public photometric data from the Kepler space mission and new radial velocity observations obtained by the SOPHIE 
spectrograph mounted on the 1.93-m telescope at the Observatoire de Haute-Provence, France. The planet KOI- 196b, with a radius of 
0.841 ± 0.032 Rji,p and a mass of 0.49 ± 0.09 Mj„p , orbits a G2V star with R^. = 0.996 ± 0.032 Rq , Mj,= 0.94 ± 0.09 Mq , [Fe/H] = 
-O.IO ± 0.16 dex, Te(f = 5660 ± 100 K and an age of 7.7 ± 3.4 Gy. KOI-196b is one the rare close-in hot-Jupiters with a radius 
smaller than Jupiter suggesting a non-inflated planet. The high precision of the Kepler photometry permits us to detect the secondary 
transit with a depth of 64*|° ppm as well as the optical phase variation. We find a geometric albedo of Ag = 0.30 ± 0.08 which is 
higher than most of the transiting hot-Jupiters with a measured Ag. Assuming no heat recirculation, we find a day-side temperature of 
Tday = 1930 ± 80 K. KOI-196b seems to be one of the rare hot-Jupiters located in the short-period hot- Jupiter desert. 

Key words, planetary systems - techniques: photometric - techniques: radial velocities - techniques: spectroscopic, star : 
individual(KOI-I96) 



1. Introduction 

Since 2007, thanks to the exoplanet-dedicated space missions 
CoRoT (Bag Un et al.|2006| and Kepler ( [Borucki et aL||2010b| l, 
the scientific community has access to very high precision 
photometry, down to a few tens of ppm with times seiies up to a 
few years. This has led to the exciting discoveries of the super- 
Earths CoRoT-lh ( |Leger et al.]|2009; Queloz et al.,.2009j and 
its twin Kepler-lOb ( BataUia et al. | 2011 1, the multi-planetary 
systems Kepler-9 ( Hoknan et al.|2010) l and Kepler-11 (JLissauerl 
et al.|20lT] l and the long-period giant planet CoRoT -9h (Deeg et 
al.|2 010')rLong space-based time series allow us to reach a very 
accurate characterization of the planetary systems, especially 
when stellar parameters can be determined by asteroseismology 
(e.g. |Batalha et al. |201 l]|Demory et al.|201 l| l. 



Moreover, high-precision space-based photometry permits 
us to find also very small effects, such as the planetary occul- 
tation (e.g. [Alonso eF al. 2009a), i.e. when the planet passes 
behind its host star, as well as the phase variation of the planet, 
i.e. the variation in brightness as the dayside of the planet rotates 
into view. The latter was found in the optical for CoRoT-lh 
([Snellen et al.|2009| l, HAT-P-7b (|Borucki et al.|2009| [Welsh et| 
al.| |2010| l, K epler-lh ([Pemoiy et al.||201 l| l and the siaper-Earth 
Kepler-\Qh ( [Batalha et al. |2011| l. 
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The first six months of Kepler data have permitted identi- 
fying 1235 planetaiy candidates around 997 stars (Borucki et 
|al.|2011 1 including twenty one fully characterized planets (with 



measured radius and mass) (e .g.|Borucki et al.|2010a[ [Latham et| 



al. 2010 [Santerne et al. 201 1 Bouchy et al. 201 1 1 and three con 



firmed planets without measured mass (e.g. Torres et al.|201 1 
Out of the remaining candidates, we have selected a few candi 
dates around stars brighter than the Kepler magnitude Kp ~ 14.7 
to be followed up with the SOPHIE spectrograph (Observatoire 
de Haute-Provence, France). One of these is the new exoplanet 
KOI-196b. In this paper we report its discovery (Sect.|2]&[3]) and 
the system characterization (Sect.|4| including an estimate of the 
planetary albedo and day-side temperature (Sect. |5]l thanks to 
the observation of the secondary and phase vaiiation (Sect. 2.3 i. 



Finally, we discuss KOI- 196b and compared it to other transiting 
planets in terms of mass, radius, period and day-side temperature 
(Sect.|6|. 

2. Kepler observations 

The Kepler Object of Interest KOI- 196 has been observed by 
Kepler since May 13, 2009 in long-cadence mode (temporal 
sampling of ~ 29.4 min). The various identifiers (ID) of this 
target, including coordinates and magnitudes are reported in 
Table [T] At the time of writing of the present paper, only 
quarters 1 & 2 (Ol & Q2) data sets are publicly available on the 
MAST archivaM They account for a total of 5 993 photometric 
measurements. 285 points were discarded by the Kepler pipeline 
( Jenkins et al.pOlOj ) as they are affected by instrumental effects. 



http://archive.stsci.edu/ 
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Table 1. KOI- 196 IDs, coordinates and magnitudes 



2.2. Primary transit modeiing 



Kepler Input Catalog (KIC) 


9410930 


Kepler Object of Interest (KOI) 


196.01 


2MASS ID 


19380317+4558539 


Right Ascension (J2000) 


19 38 03.19 


Dechnaison (J2000) 


45 58 53.76 


Kepler magnitude'"' 


14.465 ± 0.02 


„(«) 

5 


14.996 


/a) 


14.410 


;■<") 


14.238 




14.182 ±0.030 


J(b) 


13.262 + 0.022 




12.938 ± 0.024 


^(*) 


12.892 ± 0.028 


E(B - V)*"' 


0.114 ±0.100 



<"> from the Kepler Input Catalog. from 2MASS 
catalog. 



The 5708 remaining photometric measurements are displayed 
in Fig.[T] Transits with a period of ~ 1.8 days and a depth of ~ 1 
% are clearly visible in the light-curve of KOI- 196 showing no 
prominent features related to magnetic activity. 
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Fig. 1. Raw and corrected light curve of KOI- 196 as provided by 
the MAST database. 



2.1. Contamination correction 

For the following studies, we used the raw light curve (hereafter 
LC) instead of the PDC-corrected light curve as recommended 
by the Kepler Tearrj^ The flux of several nearby stars contam- 
inates the target one and has to be taken into account when 
working with the raw LC. We estimated the contamination value 
by comparing the raw and the PDC-corrected LC following the 
same methodology as described in Bouchy et al.| (201 1 ). We 
calculated it separately for each quarters. We find a contamina- 
tion value of 5.1% ± 0.06% for Ql and 3.3% + 0.04% for Q2 
which is in good agreement with the values for seasons 3 and 
0, respectively, available on the MAST database (see also Sect. 



Before performing the transit modeling, we normalized the 
transits by fitting a parabola to the ~ 1 Ih intervals of the LC 
before the ingress and after the egress of each transit in order 
to correct for any local variations. We discarded one of the 63 
available transits occuring at BJD~2455033.3 and is affected by 
one of the discontinuities seen in the LC (see Fig.[T]l. 

Transit modeling was performed following the formalism 
of Gimenez (2006i fixing the eccentricity to zero (see Sect. 
2.3 1. The seven free parameters used in the transit modeling 
were the orbital period P, the epoch of the first transit fo, the 
transit duration T14, the planet to star radii ratio k - Rp/Ri,, the 
orbital inclination /, and the two limb-darkening coefficients 
U-^ — Ua + Uh and — Ua — 

The best-fit of the primary transit was performed using the 
algorithm AMOEBA (Press et al. 1992 ) and changing the initial 
parameters with a Monte Carlo to find the global minimum of 
the;if^. We followed the procedure described in Kipping (2010|l 
to take the long cadence rate of the Kepler LC into account: the 

was computed by binning a five-times oversampled model 
LC within the long cadence rate of Kepler and compared to the 
data. The phase-folded transit and the best-fit model as well 
as the residuals between the observations and the model are 
shown in Fig. [2 We used a bootstrap procedure to estimate the 
uncertainties which consist in generating synthetic data sets 
shifting the residuals from the best-fit, adding them back to the 
best-fit model LC and fitting once again the data. We finally 
determined the 1-cr uncertainty as the 68% confidence interval 
defined as 16% above the upper (and below the lower) confi- 
dence limit of the cumulative probability. The parameters of the 
best-fit model and their 1-cr uncertainties are reported in Table[2] 
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Fig. 2. Phase-folded unbinned light curve of the KOI- 196b tran- 
sit with its best-fit model (top panel) and residuals (bottom 
panel). 



2.3. Modeling the occultation and orbital phase variations 

Space-based high-precision photometry has shown its capability 
to detect planetary occultation and phase variation for close-in 
planets (e.g. [Alonso et al.||2009a|b[ [Snellen et al.|[2009l [2010 



http://keplergo.arc.nasa.gov/PipelinePDC.shtml 



[Borucki et al.[[2009[ [Demory et ar.[[20Il l. Even if KOI- 196 
is faint, we cleaned the LC to try to detect both effects. To 
that purpose, we first removed all the transits from the raw 
LC corrected from the contamination (Sect. \2.l\ . We then 
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corrected for long-term trends separately in the Ql and the 
two first subsets of the Q2 data between BJD 2455000.0 and 
2455033.3, by fitting polynomials of the 8''' order For the 
other subsets of the Q2 LC, we just removed the offsets around 
each discontinuity. We finally de-trended the whole LC using a 
sliding median with a window of 1.5 times the orbital period, 
taking care of mirroring the LC at the beginning and the end. 
This method leads to a total dispersion of 243ppm after a 3-cr 
clipping. The obtained LC, phase-folded at the transit ephemeris 
is displayed in Fig.[3]and clearly shows an occultation as well as 
orbital phase variations. We also tested different sliding median 
windows between ~ 1.5 ■ P and ~ 2.5 ■ P and got the same result 
with a slightly larger dispersion. As expected, windows with 
extensions less than the orbital period significantly reduce the 
amplitude of the phase variations we are searching for and thus, 
our choice of L5 ■ P is the best compromise between preserving 
the signal we want to detect and minimizing the RMS of the 
residuals. 

We characterized the occultation and phase variations using 
the formalism proposed by [Snellen et al.| ( |2009 ). To model the 
occultation, we fixed the values of the orbital period, the epoch 
of transit, the transit duration, the radii ratio and the inclination 
to those derived by the transit modeling but without introducing 
the LD coefficient, as the planet is passing behind the star. We 
first binned the LC to 0.02 in phase (inset of Fig. [3]l and fitted an 
out-of- transit model with four free parameters: the phase of the 
secondary eclipse 0occ, the contrast between the planet day-side 
and stellar flux R^^y, the ratio of the night-side to day-side 



flux Fn/d and the relative stellar brightness ziev (Snellen et 
aL] 2009 | l. As for the transit modeling, we found the best-fit 
solution using the algorithm AMOEBA and changing the initial 
parameters with a Monte Carlo method to determine the global 
minimum of the x^- We also oversampled the model to take 
the bin size into account. We found the phase of the secondary 
eclipse to be ^occ = 0.499!]]™ which allow us to constrain 
ecosw = -0.002:° (|Gimenez| |2009||. Thus, we assumed 
that the eccentricity of KOI-196Bis zero. Then, we performed 
again an out-of-transit modeling to the unbinned LC, fixing the 
phase of the occultation to 0.5 and allowing the three remaining 
parameters to vary. We computed the by comparing the 
data with a five-times oversampled model binned at the Kepler 
cadence rate. Indeed, we expected that not oversampling the 
model could lead to an underestimation of the occultation depth. 
The best fit corresponding to the global minimum of the x^ is 
plotted in Fig. |3] and the derived parameters are displayed in 
Table |2] One-cr uncertainties were determined using a bootstrap 
procedure as for the primary transit (see Sect.|2.2|i. 



2.4. Centroid motion 

In order to test whether the centroid behaviour was compatible 
with the target's vicinity ( [Batalha et al.|[20T0) l, we performed 
simple 2D simulations of the photometric properties of KOI- 
196 and its five neighbours located less than 0.3' from KOI- 
196 and described in the MAST archive. All of them are fainter 
than Kp = 16.2. We used circular apertures instead of the ac- 
tual Kepler PSF; this assumption should impact only the farthest 
contaminants. It results that all configurations where one of the 
neighbours produces the transit are discarded, since this would 
imply a shift in the X and/or Y directions of 2 to 20 mpixels in 
absolute value. In our simulations, when KOI- 196 is the transit- 
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Fig. 4. Rain plots showing X (left pannel) and Y (right pannel) 
shift of the centroid during the transit as function of the normal- 
ized flux residuals. 



ing star, we still find a little centroid shift of 1-2 mpixels. The 
observed centroid positions have a mean value of in both di- 
rections, with 1-cr errors of 0.4 and 0.6 mpix in X and Y, re- 
spectively (see Fig|4]). The simple simulations of the photocen- 
ter of KOI- 196 being disturbed by its neighbours are therefore 
in agreement with the observations, within 2-cr. 



3. SOPHIE observations 

3.1. Observations and data reduction 

We performed spectroscopic follow-up observations on the 
target K OI- 196 with the SO PHIE spectrograph ( |Perruchot et al. 
2008; Bo uchy et al.|2009c I mounted on the 1.93-m telescope at 
Observatoire de Haute-Provence, France. We acquired twelve 
high-resolution spectra from March 25 to July 29, 201 1|^ using 
the High Efficiency mode (R ~ 39 000 at 550nm) of SOPHIE 
and the slow CCD read-out mode in order to minimize the 
instrumental and photon noise. Observations were done by 
keeping the SNR of the spectra constant in order to minimize 
the charge transfer inefficiency effect (Bouchy et al. 2009a i. 
Spectra were reduced with the online standard pipehne and 
radial velocities were obtained by computing the weighted 
cross-correlation function (CCF) of the spectra with a numerical 
spectral mask of a G2V star (Baranne et al. 1996; Pepe et al. 
2002). Three of the spectra were significantly affected by the 
Moon's scattered light. Their corresponding radial velocity were 
corrected using the same technique as in Santerne et al. ( |201 1\ 
and Bono mo et al.| (pOlO). Other two spectra were slightly 
affected by the Moon scattered light and were not corrected 
since the Moon light affected the measurement at the level of 
a few ms \ which is compatible with the noise added by the 
correction. 



These radial velocities are listed in Table [3] and displayed 
in Fig.|5]&[6] They show a clear variation compatible with the 
reflex motion of the parent star KOI- 196 due to a planetary 
companion in phase with the Kepler ephemeris. Since we found 
a secondary eclipse at phase ~ 0.5 (c.f. Sect. 2.3 i and that 
circularization time-scales for such close-in planets are very 
short ( Mardhng|2007 1, we assumed the orbit to be circular and 
fitted a corresponding keplerian to the data. We found a best fit 
with a semi-amplitude = 85 + 1 1 m s"' and a ctq-c = 24 m s"' 
which is comparable with the mean radial velocity uncertainty 



prog. ID: IIA.PNP.MOUT 
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Fig. 3. Unbinned phase-folded light curve of KOI- 196 with the best-fit model. The inset displays a zoom of the out-of-transit LC 
binned to 0.02 in phase with the best-fit model of the planet phase variation and occultation (inset, top panel) and its residuals (inset, 
bottom panel). Out-of-transit residuals binned to 0.02 in phase have a dispersion of 22 ppm. 



{cr„) = 20 ms 

From the SOPHIE CCF parameters, one can estimate the 
V sin and [Fe/H] using equations described in Boisse et al. 



( |2010| l. Using the FWHM and contrast of the CCF and assum- 
ing {B — V) - 0.63, we find a i/sin/* = 4.1 ± 1.0 kms ' and 
[Fe/H] = -0.17 ± 0.15 dex. 




0.4 0.6 
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Fig. 5. Phase-folded radial velocity curve with our best-fit (top 
panel) and residuals (bottom panel). 
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Fig. 6. Radial velocity curve with our best-fit (top panel) and 
residuals (bottom panel). 



and radial velocities or different RV amplitude when processing 
the CCF with diff'erent spectral type mask templates ( ]Bouchy| 
|et al. ,2009b) . We first fitted the different RV datasets processed 
with K5 and FO masks with a circular orbit model and found 
that /Tks = 84 + 11 ms"', cro_c = 20 ms"' and K^q = 85 ± 14 
ms"\ (To-c = 32ms"' which is in good agreement with the 
amplitude derived with the G2 mask. 



3.2. Blend analysis 

If the RV signal is mimicked by a diluted blended binary one 
can expect to see some correlation between the bisector span 



Bisector span is also a key CCF diagnostics which is very 
sensitive to stellar activity and additional blended stellar compo- 
nent ( Queloz et al.j2001 ; Bouchy et al. 200 9b, Santerne et al., in 
prep.). In order to assess the possibility that the radial velocity 
variations are not caused by a blended binary, we measured 
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the bisector spans which are listed in Table [3] and plotted in 
Fig. |7] They do not reveal any significant variation within 1-cr. 
Moreover, neither radial velocities, nor their residuals present a 
significant correlation with the bisector These two checks allow 
us to secure the planetary nature of KOI- 196b. 
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Fig. 7. CCF bisector as function of radial velocity (left panel) 
and radial velocity residuals (right panel). 



an upper value of the reddening provided in the MAST database. 

Finally, we used the stellar density derived by the transit 
modeling combined with the Tgff from the spectral analy- 
sis to determine the star's fundamental parameters in the 

1 /3 

(Teff , M^' //?*) space. The position of the star in the H.- 
R. diagram was compare d to STAREVOL evolution tracks 



( Turck-Chieze et aL]2010| l minimizing the as described in 



Santerne et al. (2011 1. We found a main-sequence solution 
with = 0.94 ± 0.09 Mq and = 0.966 + 0.032 R0 . 
Main-sequence solutions are consistent with the absence of 
Lithium in the spectrum of the star and the low v sin /* value 
we derived. We note that the inferred gravity surface log 
g= 4.44 ± 0.10 is in good agreement with the spectroscopic 
value within 1-cr as well as the one published by Borucki et 
[aLl ( [20TT) : log g=4.5. This solution led to an age of 7.7 + 3.4 Gyr. 

These stellar parameters yield limb darkening coefficients in 
the Kepler bandpass of m+ - 0.66 + 0.03 and m_ = 0.15 + 0.03 
(using interpolated values from Sing 2010). They are in good 
agreement with those obtained from the transit modeling within 
2-cr and 1-cr for m+ and m_, respectively. 



4. System parameters 

To perform the spectral analysis of the KOI- 196b host star, 
we used the four SOPHIE spectra took before June 2011 
that were not affected by Moon reflected light at the time of 
the observations. Since such individual spectra have too low 
signal-to-noise ratio to allow a proper analysis, we co-added 
them after correcting of the barycentric Earth radial velocity 
projected through the line of sight and of the star's one. The 
resulting co-added spectrum has a signal-to-noise of 46 per 
element of resolution at thecontinuum at 5550 A with a spectral 
resolution of ~ 39000. We used the semi-automatic software 
package VWA ( |Bruntt et aL][2002l [2008) to derive the st ellar 
photospheric parameters as described in |Bruntt et al. (20101. 



et al. 



Using the same methodology as presented in Bruntt 
( 2010j ) we derived the projected rotational velocity 
i/sin;* = 4.5 ±1.5 kms"' from the analysis of a few 
isolated spectral lines. This value is in good agreement with the 
V sin ii, estimated from the CCF assumed that the true value of 
the V sin could be even lower since we are here limited by the 
spectral resolution of SOPHIE. For the atmospheric parameters 
we found Teff = 5660 ± 100 K, [Fe/H] = -0.09 + 0.16 dex 
and log g= 4.35 ± 0.2. It is known that the log g is the most 
uncertain parameter derived from spectral analysis. The quite 
low signal-to-noise of our co-added spectrum combined with 
the moderate spectral resolution do not permit us to estimate 
this parameter more accurately. The derived stellar parameters 
are reported in Table |2] They are in a good agreement within the 
error bars, with the estimates published in Borucki et al. ( 201 l| l 
Teff - 5585 K and log §=4.51 with uncertainties up to 30%. 

From infrared photometry J and K available in the 2MASS 
archive, we computed the photometric Teff = 6000 + 150 K fol- 
lowing ^asagrande^eF^T] ( 2010|, assuming [Fe/H] of -0.10 dex 
and taking a reddening of E(J-K)=0.059 ( jCardelli et a"L|[T989] l 
into account. This photometric estimation is at 2-cr from the 
one determined by spectral analysis and might be explained by 



From the adopted stellar fundamental parameters, com- 
bined with results from the transit modeling and the anal- 
ysis of the radial velocity observations, we find a mass 
and a radius of the planet of Mp - 0.49 ± 0.09 Mjup and 
Rp - 0.841 + 0.032 Rjup . The inferred mean density of the 
planet is pp = 1.02 + 0.18 g.cnT^ . 



5. Geometric albedo characterization 

The high accuracy of the Kepler photometry as well as the large 
number of orbital periods covered during Ql & Q2 allowed us 
to clearly detect both the secondary transit and the phase varia- 
tion in the optical (a 5-cr detection, see Sect. \23\ . Assuming a 
pure reflecting planet without thermal emission we found a geo- 



metric albedo of Ag - 0.33!|5|5g in the Kepler bandpass (Rowe et 
[aLl|2006| l. Even in the optical, part of the observed light from 
the planet could be due to thermal emission (Lopez-Morales 
'& Seager"2007J [Snellen et aL]|2009] l. The observed occultation 
depth is thus a combination of the reflected light and the thermal 
emission: 



^occ — 



Plh 



^ Pref 



(1) 



where is the thermal flux generated by the planet and Fp^.^ 



is the reflected light. Using eq. 14 in Rowe et al. (2006i the re 



fleeted component of the occultation depth expressed as: 



(2) 



for which Rp/a is well constrained by the modeling of the 
primary transit. 

To estimate the fraction of thermal emission that contributes 
to the planetary occultation, we computed the ratio between the 
flux of the star and the expected thermal flux of the planet as- 
suming a blackbody emission. If the planet is tidally locked, its 
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Table 2. Star and planet parameters. 



Ephemeris 


Planet orbital period P [days] 


1.855558 ±7.10-" 


Transit epoch T,, [BJD - 2454900] 


70.1803 ±0.0003 


Results from radial velocity observations 




Orbital eccentricity e 


(fixed) 


Semi- amplitude [ m s"' ] 


85 ± 11 


Systemic velocity V,. [kms ] 


-27.066 ± 0.007 


O-C residuals [ms"'] 


24 


Fitted transit parameters 




Radius ratio k = Rp/R^ 


0.0895 ± 0.0019 


Orbital inclination / [deg] 


88.3 ± 0.7 


LiailML QUlaLlUll 1 [4 j^llj 




limb darkening coefficient u+ 


0.98± 0.17 


limb darkening coefficient u_ 


0.16 ± 0.20 


Deduced transit parameters 




Scaled semi-major axis a/R^, 


6.43 ± 0.05 


Impact parameter b 


0.19 ±0.07 


A^i'V^* [solar units] 


1.012 ±0.009 


Stellar density [g cm^^] 


1.46 ±0.04 


Fitted and deduced Out-of-Transit parameters 




phase of secondary transit (p^cc 


"•^'^-0.011 


planet day-side to stellar flux ratio R^-^y [ppm] 


64:;° 


planet night-side to day-side flux ratio F^jq 




relative stellar brightness z\c^ 


0.99993 ± 0.00001 


e cos a) 


-0 002+°°°'* 


Spectroscopic parameters 




Effective temperature Teff [K] 


5660 ± 100 


MetalUcity [Fe/H] [dex] 


-0.09 ±0.16 


Stellar rotational velocity u sin [ km s"' ] 


4.5 ± 1.5 


Spectral type 


G2V 


Stellar physical parameters from combined analysis 


Star mass [Mq ] 


0.94 ± 0.09 


Star radius R^, [Rg ] 


0.966 ± 0.032 


Surface gravity log g" [dex] 


4.44 ±0.10 


Age of the star t [Gy] 


7.7 ± 3.4 


Distance of the system [pc] 


730 ± 70 


Planetary physical parameters from combined analysis 


Orbital semi-major axis fl* [AU] 


0.029 ± 0.001 


Planet mass M,, [Mjup ] 


0.49 ± 0.09 


Planet radius i?,,[Rjup ] 


0.841 ± 0.032 


Planet density p,, [g cm^^] 


1.02 ±0.18 


Geometric albedo A„ 


0.30 ± 0.08 


Day-side temperature Tjay [K] 


1930 ±80 



" derived from and R*. 
* derived from a/R* and R*. 
considering no atmospheric thermal circulation. 



day-side temperature can be calculated using the following ex- 
pression ( Cowan & Agol|2011[ l: 



"12^ 



(3) 



where the Bond albedo is fixed to zero assuming black- 
body emission and e is the circulation efficiency that is allowed 
to vary from e = 1 for a full redistribution of the energy budget 



from the day-side towards the night-side to e 
circulation in the atmosphere. 



for no heat 



We used the 



-CK 



Castelli & Kurucz 



(2004 1 atmospheric mode|^ 
, of a solar-like star to estimate the stellar flux and integratec 
over the Kepler bandpass Q,j^ Then, we integrated the Planck 
function at the day-side brightness temperature of the planet 
over the Kepler bandpass. The occupation depth due to thermal 
emission is the ratio between the two fluxes, multiplied by the 
surface areas of the planet and the star, respectively : 



Pth 



(4) 



where h is the Planck constant, ks, the Boltzmann constant and 
c the speed of light in the vacuum. 

Fig.[8]displays the geometric albedo in the Kepler bandpass 
as a function of the brightness temperature of the planet as 
derived from the previous equations. We found that the thermal 
emission of the planet is negligible within the Kepler bandpass 
and we found a geometric albedo of Ag - 0.3 ± 0.08. Using 
this range of values for the albedo and assuming that the geo- 
metric albedo related to the Bond albedo through the relation: 
Ab = 2/3Ag for a perfect Lambertian sphere (e.g. Rowe et al 



20061, and assuming that the energy redistribution is not e 
cient (s = 0) for such a strongly irradiated planet, we estimated 
the day-side temperature of the planet of Tday - 1930 + 80K. 




1800 2000 2200 2400 

Day-side brightness temperature [K] 



Fig. 8. Geometric albedo as function of planet day-side bright- 
ness temperature. Dashed lines represent the brightness tempera- 
ture of the planet assuming zero Bond albedo and a perfect heat 
redistribution , e = 1 ; T=1590 ± 60 K (left) and no thermal 
circulation in the atmosphere e = ; T=2030 + 60 K (right). 
The grey band covers the albedo values allowed by the 1-cr un- 
certainty on the occupation depth which is assumed to be the 
dominant uncertainty. 



http://www.stsci.edu/hst/observatory/cdbs/castelli_kurucz_atlas.html 
http://keplergo.arc.nasa.gov/CalibrationResponse.shtml 
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6. Discussion 



6.1. A non-inflated low-mass hot-Jupiter 

With an orbital period of 1.855558 ± 7.10"^ d, a radius of 0.841 



0.032 Rjup and a mass of 0.49 ± 0.09 M. 



Jup 



KOI-196b is 



one of the rare objects with a very short orbital period (P < 2 
d) and radius and mass smaller than Jupiter Fig. [9] shows 
KOI- 196 compared with other planets in terms of orbital period, 
planetary radius and mass. Only HD212301b ( |Lo Curto et 
[al7 2006 ) exhibits similar characteristics in terms of mass and 
orbital period. However as the orbital inclination of this planet 
is still unknown, the true mass of HD212301b could be likely 
much higher. Using tables in Fortney et al. ( 20071, we estimated 



that the core mass of KOI- 196b possesses between 50 M® 
and 100 of heavy elements. This means that between ~ 
1/3 and ~ 2/3 of KOI-196b's mass is concentrated in its core. 
Another high-density short-period planet with radius smaller 
than Jupiter's one is HD 149026b ( |Sato et al.|[2005l ). Out of the 
141 transiting planets discovered so far, only KOI- 196b and 
HD 149026b are located in the hot- Jupiter desert with orbital 
period shorter than 3 days and with radius or mass smaller than 
Jupiter. This desert could not be explained by an observational 
bias since their short orbital periods make their detection easier 
These two objects might share the same formation and/or 
evolution process ( |Broeg & Wuchterl|2007 1. 



Fig. 10 displays the planetary radius as function of the 



maximum day-side temperature assuming no heat circulation 
(i.e. r^^" = Teif / y/a/Ri,). There is a clear correlation between 
the radius and the expected maximum day-side temperature for 
hot-Jupiters less massive than 2 M jup as expected theoretically 
( Burrows et al.| 2007 Liu et al. 2008| l. Once again, in the 
Saturn- Jupiter domain, KOI- 196b and HD 149026b appear to 
be outliers. This might be explained by an unusual high Bond 
albedo in comparison with other transiting planets as it is 
observed for KOI- 196b. Indeed, it is expected that planets with 
a higher albedo are cooler (see eq. [3]i and thus, are less inflated 
(seeFig.fTOli. 



E 1.0 



■ non transiting planets 
• transiting planets 
^ KOI-196b 



■ ♦ 



* ■ ■ 



• transiting planets 

• KOI-196b 

• HD149026 b 



4 I I » HD14; 



3 4 
orbital period [d] 



Fig. 9. Planetary mass (upper panel) and radius (lower pannel) 
as a function of orbital period for exoplanets with orbital period 
less than 6 days (source: Wright et al. 20111. Transiting exo- 
planets are the red circles and RV planets are the blue squares. 
We note that for non-transiting planets, the minimum M sin / is 
represented here. KOI-196b is marked by a black diamond and 
HD149026b is marked with a magenta hexagon. 



High mass-loss can explain this unusual low mass and low 
radius for a such short-period hot- Jupiter (Ehrenreich & Desert 
201 l) l. Using eq. 15 and 22 in jLecaveher Des Etangs, (2007( ) 
and assuming a efficient extraction of the gas by the extreme 
ultraviolet incident flux rj - 1, we estimate a current escape rate 
for KOI- 196b of ~ 0.0005 Mj^p / Gyr Thus, mass-loss is a very 
unlikely process to explain the low mass and small radius of 
KOI- 196b. 

Further studies of internal structure, taking the high albedo of 
KOI- 196b into account, are required to explain the small radius 
for this close-in planet. 

6.2. A high albedo for a hot-Jupiter 

KOI- 196b is the third hot Jupiter, discovered so far with a such 
high albedo. The first ones were Kepler-lh (Latham et a l. [2010)1 
wit h Ag = 0.32 + 0.03 (Kipping & Bakos 2011; Demory et 
[aL][20Tl] l and HA T-P-7b IjPk et al. 2008 ) with an albedo of 
Ag = 0.58 + 0.05 ( [Welsh et al.|2010| l. Other giant planets with 
an observed occultation have a very low albedo in the optical 
Ab < 0.1 ( Cowan & Agol 2011 ). Theoretical studies previously 
predicted the low-albedo of close-in planets (e.g. Sudarsky et 
[ar]|2000| l due to the high opacity of Na and K at optical wave 



lengths, [bemory et al.| ( |20lT] l suggested that a depletion of such 
elements in the high atmosphere of Kepler- 7b compared to solar 
abundances can explain its high albedo. Sudarsky et al. ( 2000[ l 
have shown that a high albedo can also be explained by the pres- 
ence of silicate layer high in the atmosphere. Finally, synthetic 
spectra of the thermal emission of hot-Jupiters from [Fortney et] 
|aLl ( [2008) suggests that thermal emission could be significantly 
different from a blackbody emission in the optical. In that case, 
contribution of the thermal light in the Kepler bandpass can be 
significantly higher, leading to a lower value of the geomet- 
ric albedo of KOI- 196b. Infrared observations of the secondary 
eclipse and, if possible, the thermal phase variation of KOI- 196b 
are needed to constrain the day-side and night-side temperature 
as well as to constrain unambiguously its geometric albedo and 
heat circulation ( [Cowan & AgoipOTT i. 



7. Conclusion 

Thanks to publicly-available Kepler photometry and new 
high-resolution SOPHIE spectroscopic observations, we es- 
tablished the planetary nature of the Kepler Object of Interest 
KOI-196.01, now caUed "KOI- 196b", one of th e 1235 Kepler's 
candidates published by Borucki et al. ( 201 l| l. KOI-196b is a 
0.49 + 0.09 Mjup - 0.841 + 0.032 Rjup hot-Jupiter orbiting a 
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Table 3. SOPHIE measurements of KOI-196. 



4 KOI-196 b 
• HD149026 b 




BJD 


RV 


±lcrrv 


v 

' spun 


S/N/pix 


-2 454 900) 


[kms-'] 


[kms-'] 


[kms-'] 


@550nm 


746.63646 


-27.063 


0.017 


-0.013 


16.8 


749.65408 


-27.141 


0.016 


0.009 


17.0 


750.65676 


-27.024 


0.025 


-0.007 


14.7 


760.60669 


-27.110 


0.018 


-0.035 


17.5 


761.61891 


-27.005 


0.015 


0.026 


19.1 


768.62681 


-27.130 


0.027 


-0.033 


14.3 


769.63221 


-26.998 


0.018 


-0.022 


19.3 


770.62035 


-27.053 


0.021 


0.022 


16.6 


854.48407 


-27.005 


0.024 


-0.015 


13.7 


868.45141 


-27.169 


0.019 


-0.051 


16.8 


871.45738 


-26.957 


0.025 


0.021 


17.4 


872.59440 


-27.082 


0.020 


-0.009 


17.9 



1000 1500 2000 2500 3000 
Maximum day-side temperature [K] 



Fig. 10. Planetary radius vs estimated maximum day-side tem- 
perature assuming no heat circulation for all transiting planet 
(source: Wright et al.|2011 1. Transiting planets with mass Mp ^ 



port from the French National Research Agency (ANR-08-JCJC-0102-01). A.S. 
is grateful to Cilia Damiani for all the fruitful discussions. A.S.B. is supported 
by CNES. 



2 Mjup are the blue squares. Transiting planets with mass in be- References 
tween Mp = 0.1 Mj^p and Mp - 2 Mjup are the red circles. 



Transiting hot-Neptunes and Super-Earths are the green trian- 
gles. KOI-196b is marked by a black diamond and HD149026b 
is marked with a magenta hexagon. 

slightly metal-poor G2V star in a 1.855558 + 0.000007d period. 
This planet is one of the rare close-in hot-Jupiters with a radius 
and mass lower than Jupiter suggesting that it is a non-inflated 
planet. We note that only 6 Kepler candidates with priority 
2, including KOI- 196b, are expected to be hot-Jupiters with 
an orbital period shorter than 3 days and radius in between 
0.5 Rjup and 0.9 Rjup : KOI-102.01, KOI-183.01, KOI-356.01, 
KOI-801.01 and KOI-883.01. Focusing follow-up efi^orts on 
them can help to specify the boundary of the short-period low 
mass and small radius hot- Jupiter desert, if exists. 

Using the high-quality and long time-series of the Kepler 
data, we detected the occultation as well as the optical phase 
variation of KOI- 196b, orbiting a quite faint star (my ~14.6). 
From the occultation depth and assuming a blackbody planetary 
thermal emission in the Kepler bandpass, we estimated a 
geometric albedo of the planet to be 0.30 + 0.08, which is 
markedly higher than the geometric albedos observed for most 
of the known hot-Jupiters ( Cowan & Agol|[201 1[ ). This leads 



to a day-side temperature of Tday - 1930 + 80 K assuming no 
thermal redistribution in the atmosphere. Lower values for the 
albedo are still possible but would indicate a significant contri- 
bution of the planetary thermal emission in the optical. This can 
be unambiguously confirmed with infrared observations of the 
planetary occultation. 

This detection demonstrates once again the efficiency of 
SOPHIE, a dedicated instrumentation on a 2-m class telescope 
for the ground-based follow-up of space mission like Kepler, 
CoRoT or the ESA M-class mission PLATO, if selected. 
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